Tachyon cosmology, supernovae data and the Big Brake singularity 
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We compare the existing observational data on type la Supernovae with the evolutions of the 
universe predicted by a one-parameter family of tachyon models which we have introduced recently 
in paper [7|. Among the set of the trajectories of the model which are compatible with the data there 
is a consistent subset for which the universe ends up in a new type of soft cosmological singularity 
dubbed Big Brake. This opens up yet another scenario for the future history of the universe besides 
the one predicted by the standard ACDM model. 

PACS numbers: 98.80.Cq, 98.80.Jk, 98.80.Es, 95.36.+X 
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I. INTRODUCTION 

The discovery of cosmic acceleration [l| has stimu- 
lated the study of different models of dark energy Q 
which may be responsible for such a phenomenon. Mod- 
els of dark energy include those based on different per- 
fect fluids, having negative pressure, on minimally and 
non-minimally coupled scalar fields and on fields having 
non-standard kinetic terms [!, H| . The latter ones include 
as a subclass the models based on different forms of the 
Born-Infeld-type action, which is often associated with 
the tachyons arising in the context of string theory ||. 
Tachyonic models with relatively simple potentials were 
confronted with observational data in [6(. Compared to 
the standard Klein-Gordon scalar field cosmological mod- 
els the dynamics of tachyon models can be much richer 
due to the non-linearity of the dependence of the tachyon 
Lagrangians on the kinetic term of the tachyon field. 

In a recent paper Q a particular one-parameter fam- 
ily of tachyon models was considered, which has revealed 
some unexpected features. At some values of the param- 
eter of the model a long period of accelerated quasi-de 
Sitter expansion is followed by a period of cosmic decel- 
eration culminating, after a finite time, in an encounter 
with a cosmological singularity of a new type, which was 
named Big Brake. This singularity is characterized by 
an infinite negative value of the second time derivative 
of the cosmological radius of the universe, while its first 
time derivative and the Hubble variable vanish, and the 
radius itself acquires a finite value. This singularity be- 



longs to the class of soft (sudden) cosmological singular- 
ities @, H, [n| which have been rather intensively studied 
during the last years. Here it is worth mentioning that 
in the context of the scrutiny of candidates for the role of 
dark energy, some other singularities attract the atten- 
tion of cosmologists. Among them a special place occu- 
pies the Big Rip singularity BjJ, arising in some models 
where phantom dark energy [12| is present. The possibil- 
ity of existence of a phase of contraction of the universe, 
ending up in the standard Big Crunch cosmological singu- 
larity was also considered in the literature [r| ■ Recently, 
w-singularities were also proposed (l4| . 

An attractive peculiarity of the tachyon model studied 
in paper 0] is the fact that there the Big Brake singu- 
larity is not put in "by hands" , but arises naturally as a 
result of the cosmological evolution, provided some ini- 
tial conditions are chosen. Therefore it is a consequence 
of the dynamics, rather than a pure kinematical possi- 
bility. Such evolution leading to the Big Brake coexists 
with another type of evolution describing an infinite ex- 
pansion of the universe. In other words, a small change 
of initial conditions can have drastic consequences for 
the future of the universe. Actually, in spite of it being 
somewhat exotic, we show that the cosmological model 
does not contradict observations. To this aim we com- 
pare the cosmological evolutions predicted in Q with the 
data coming from the supernovae type la observations. 
We select the compatible initial conditions by studying 
the backward evolution in comparison with the luminos- 
ity - redshift diagrams for the supernovae type la stan- 
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FIG. 1: (Color online) Phase portrait evolution for k > 
(k = 0.44). 



dard(izable) candles. Then, choosing initial conditions 
which are compatible at the la level with the data, we 
study the forward evolution and show that a decelera- 
tion period following the present accelerated expansion 
is possible, and when it is so, we estimate how long it is 
expected to last. 

The structure of the paper is the following. In Sec. II 
we introduce the model and its basic equations; in Sec. 

III we find a subset of initial conditions which are com- 
patible with the observational data by integrating numer- 
ically the dynamical equations backwards in time; in Sec. 

IV we study numerically the cosmological evolutions for 
the selected initial conditions by numerical integration 
forward in time. We end with some concluding remarks. 



II. TACHYON COSMOLOGICAL MODEL 

We consider the flat Friedmann universe with the met- 
ric ds 2 = dt 2 — a 2 (t)dl 2 , filled with a spatially homo- 
geneous tachyon field T evolving according to the La- 
grangian 



L = -V(T)dl-g Q0 T 2 . 



(1) 



The energy density and the pressure of this field are, 
respectively 



V{T) 

V\-f 2 



and 



p = -V{T)V1-T 2 . 

The equation of motion for the tachyon is 

f at V T 
+3— + -£■ = 0. 

1 - T 2 a V 



(2) 



(3) 



(4) 



We consider the following tachyon potential V(T) @: 

A 



V(T) 



sin 2 (ly/A(l + k)T) 



'l-(l + fc)cos 2 (^y/A(l + k)Tj, 



(5) 



where A is a positive constant and — 1 < k < 1. 

Taking into account the Friedmann equation H 2 = e, 
where the Hubble variable H is defined as H = a/ a, and 
the Newtonian constant is normalized as 8ttG/3 — 1, we 
obtain the following dynamical system: 



T = s, 



(0) 



s = -3^(l- S 2 ) 3 / 4 S -(l- S 2 )-^. (7) 

When the parameter k is negative, the evolution of the 
system @-([7]) is confined inside the rectangle 



1< s < 1, 



< T < 



2tt 



3^A(l + fc) 
The system has only one critical point: 



3^A(l + fc) 



, so = 0, 



(8) 



(9) 



(10) 



which is an attractive node corresponding to a de Sitter 
expansion with Hubble parameter 



Ha 



A. 



(11) 



All cosmological histories begin at the Big Bang type 
cosmological singularity located on the upper (s = 1) or 
lower (s = —1) side of the rectangle (|!Sj)— (J9j> , the individ- 
ual history being parametrized by the initial value of T 
satisfying the inequality © ■ They all end up in the node 

m- 

In the case k > the situation is more complicated. 
First of all, the real potential V is well-defined only in 
the interval 



where 



T 3 <T< T 4 , 



2 1 
la = — larccos 

3^/(1 + Jfc)A VT+k 



Ta 



3^/(1 + fc)A 



7r — arccos- 



(12) 



(13) 



(14) 
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FIG. 2: (Color online) The fit of the luminosity distance vs. redshift for k — —0.4 (upper left), —0.2 (upper right), (middle 
left), 0.2 (middle right), 0.4 (lower left), 0.6 (lower right), in the parameter plane (yo, wo — 1/ (l + So)). The white areas 
represent regions where the bounds on the model are not satisfied. The contours refer to the 68.3% (1<t) and 95.4% (2a) 
confidence levels. For increasing values of |fc| < 1 the well- fitting regions are increasingly smaller. The colour code for \ 2 is 
indicated on the vertical stripes. 



The dynamical system ©-([T]) has three fixed points: 
the node ((TO)) and the two saddle points with coordinates 



Ti = 



3v/(l + fc)A 



1 - k 

arccos-W -, Si = 0, 

1 + k 



(15) 



and, respectively, 
T, 



n aiccost / % I , S2 — 0, (16) 



3V(1 + *)A V'' "Vl + fc, 

which give rise to an unstable de Sitter regime with Hub- 
ble parameter Hi = yj ^^yj^ > Ho- 
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The most striking feature of the model under consid- 
eration with k > consists in the fact that now the 
cosmological trajectories do cross the corners of the rect- 
angle (f8"| . (fT2"|) . Indeed, the direct analysis of the sys- 
tem of differential equations in the vicinity of the points 
P, Q, Q' and P' (see Fig. [1]) shows that these points are 
not singular points of the system Moreover, there 
is no cosmological singularity in these points 0]. That 
means that the cosmological evolutions must be contin- 
ued through them. An apparent obstacle to such a con- 
tinuation is the fact that the expression under the square 
root in the formula for the potential ([5]) changes sign 
when T becomes smaller than T3 or greater than T4. 
However, the expression under the square root for the 
kinetic term yl — s 2 also changes sign at the same time. 
Then, since the Lagrangian of the theory is the product 
of these square roots, these simultaneous changes of sign 
leave the Lagrangian and the corresponding expressions 
for the energy density @ and the pressure real. The 
equation of motion for the tachyon field (0} also conserves 
its form. The sign, which we prescribe for the product (or 
for the ratio) of the square roots is uniquely determined 
by the Friedmann equation. In analyzing the behavior of 
our dynamical system in the regions where \s\ > 1 it is 
convenient to use the new potential 



W(T) 



A 



v/A(i + fc)r) 



'(i 



fc)cos 2 ( -y/A(l + k)T 



1. 



(17) 



and to substitute in all expressions the term 1 — s 2 by 
s — 1. In doing so the energy density and pressure have 
the form 

W{T) 



e = 



1 



and 



P : 



W(T)^s 2 - 1, 



(18) 



(19) 



being both positive. 

The procedure of continuation of the trajectories 
through the corners of the rectangle is described in de- 
tail in Q. Here, for the convenience of the reader we 
reproduce the phase portrait of the dynamical system 
from [7J with some brief comments. The rectangle in 
the phase space (T, s) should be complemented by four 
infinite stripes (see Fig. Q]). The left upper stripe (the 
right lower stripe) corresponds to the initial stages of the 
cosmological evolution, while the right upper stripe (the 
left lower stripe) corresponds to the final stages. There 
are five classes of qualitatively different cosmological tra- 
jectories. The trajectories belonging to classes I and II 
end their evolution with an infinite de Sitter expansion, 
while the trajectories of classes III, IV and V encounter 
a Big Brake singularity. The curves u, £, r, -0 and \ are 
separatrices, dividing different classes of trajectories. 



We end this section with the following remark. Like the 
other tachyon or DBI cosmological models (for example, 
models displaying the power-law or exponential poten- 
tials) the model based on potential ([5]) possesses a wide 
class of cosmological evolutions ending up in an infinite 
accelerated expansion. In addition, for small values of T, 
this potential behaves as 1/T 2 , a behavior which has been 
widely studied in the literature. So far, so good. On the 
other hand, because of the more complicated structure 
of the potential ([5]) , our model exhibits another class of 
trajectories with a qualitatively very different behavior 
and, in our opinion, this is precisely the feature which 
makes it particularly interesting. 



III. THE TACHYON COSMOLOGICAL MODEL 
AND COMPARISON WITH SUPERNOVAE 
TYPE IA OBSERVATIONAL DATA 

In this section we select, at the confidence level of lcr, 
and for a given choice of values of the parameter k the 
set of initial conditions (z = 0) for the system ©-((7]), 
which are compatible with the supernovae type la data 
taken from paper [l5j . To this purpose, for the numerical 
analysis of the model it is convenient to rescale the rel- 
evant variables introducing the following dimensionless 
quantities: 



where Hq is the present value of the Hubble parameter 
Ho = H(z = 0). In addition we find it convenient to 
replace the variable T with the new variable 



cos[-^n A (l + k)T 



also to switch from the time derivative to the derivative 
with respect to the redshift z: 



(21) 



and denote d/dz with a prime. 

Then, the system of equations ©-(O in terms of the 
new variables H , s, y (all depending on z) becomes: 



H z = 



V 



(l + z)Hs' = 



(1- S 2 ) 1/2 ' 

2y'(l + z)H 
3y/ n A (1 + k) (1 - y 2 ) 

+ (l-s 2 )^, 



(22) 
(23) 

(24) 
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where V and Vt are given by 
n A [1 - (1 + k) y 2 



V = 



i-y 2 



3fi Av / n A (l + k)y [k - 1 + (1 + fc) y 2 



2 (l_ y 2)3/2 [1 - (1 + fc)y 2 ] 



1/2 



(25) 



.(26) 



ables s and ?/ satisfy the constraint 



«(0)=± 



\ 



1 - 



"a 



(l + fc)y(or 



[l-y 2 (0)r 



We can avoid double coverage of the parameter space (the 
model being invariant under the simultaneous change of 
signs yo — » — yo and Sq — ► — So) by replacing so by the 
new variable 



1 



w 



1 



(27) 



The luminosity distance function for a flat Friedmann 
universe 



d L (z) = (1 + z) 



dz* 



o H(z*) 



(28) 



gives for the dimensionless luminosity distance = 
Hod^ the equation 



d L 



1 



(29) 



We are now in a position to compare our model with 
the Supernovae type la data [HI . 

Following Ref. [16[ we introduce the distance modulus 
type quantity 5 log 10 dL (z)+M, with M a constant offset 
between the data and the theoretical expression. The 
comparison involves computing 



A' 



X 2 = X! ~2 5 lo Sio d T P ( z ») -M-5 log 10 d L (zi) 

(30) 

where the sum is over the supernovae in the data set and 
<ii are the experimental errors in 5 log 10 (i L xp (zi). The 
distance luminosity function d^ (z) depends on the initial 
condition yo — 2/(0) an d sq = s (0). We minimize this 
expression with respect to M obtaining 



L 

= ~~D 



(31) 



with 



N 



L 



D 



^ 2 
N 1 



5 logio dT P {z^ - 5 log 10 d L (z^ ,(32) 



(33) 



In TableUare listed the values yj, j = 1, 2, 3, 4, of the 
variable y corresponding to the values Tj of the variable T 
given in formulas p3 ]) - (fT6|l for the chosen positive values 
of k. 



TABLE I: The values of yj (corresponding to the Tj) for some 
positive values of k. 



vari- 


k 


0.2 


0.4 


0.6 




2/1,2 


±0.816 


±0.655 


±0.500 




J/3,4 


±0.913 


±0.845 


±0.791 



Since the expansion of the present day universe is ac- 
celerated the pressure is negative, hence |so| < 1. There- 
fore, the initial point in the phase diagram (T, s) should 
lie inside the rectangle (T 3 < T < T 4 , \s\ < 1), (see Fig. 
[IJ. Thus the bounds on the model are not satisfied in 
the ranges y < y A and y > 2/3- 

In Fig [2] we represent the values of x 2 in the parameter 
plane of the initial conditions (t/o = 2/(0), wq = w(0)), for 
the choices k — 0, ±0.2, ± 0.4 and 0.6. The contours 
represent the 68.3 (la) and 95.4 (2a) confidence levels 
and the white areas are unallowed regions. 



IV. FUTURE COSMOLOGICAL EVOLUTIONS 

In this section, in order to investigate the possible fu- 
tures of the universe within the tachyon cosmological 
model, we evolve numerically the model forward in time 
starting from the parameter range (wq, yo) of initial con- 
ditions for which the fitting with the supernovae data is 
within la (68.3%) confidence level. We do this by nu- 
merical integration of equations of motion from z = 
towards negative values of z. 

The results of these computations, corresponding to 
the six values of k chosen earlier, are displayed in Fig. [3] 
in the space (w = (1 + s ) , y, z). The evolution curves 
start from the allowed region (wq, yo) in the plane z — 0. 
The final de Sitter state is characterized by the point 
(wds = 1- yds — 0, Zds = — 1), the Big Brake final state 
by points (w B b = 0, -1 < ysB < 0, -1 < z B b < 0). 

Whereas all trajectories with k < end up eventually 
into the de Sitter state, those with k > can either evolve 
into the de Sitter state or into the Big Brake state, de- 
pending on the particular initial condition (wo,yo)- The 
fraction of curves eventually meeting a Big Brake in- 
creases with increasing k. This is clearly seen in Fig. 
[3] from the relative sizes of the la subdomains belonging 
to these two regimes, which are separated by a line. 

For all future evolutions encountering a Big Brake 
singularity we have computed the actual time tsB it 
will take to reach the singularity, measured from the 
present moment z = 0, using the equation (Hot)' = 
-H- 1 (1 ± Z y l ). The results are shown in Tables UllV] 
In the tables the parameter values at which the pressure 
turns from negative to positive are also displayed. 



»=i 
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TABLE II: Properties of the tachyonic universes with k — 0.2 
which (a) are within la confidence level fit with the type la 
supernova data and (b) evolve into a Big Brake singularity. 
Columns (1) and (2) represent a grid of values of the allowed 
model parameters. Columns (3) and (4): the redshift z* and 
time t, at the future tachyonic crossing (when s = 1 and 
the pressure becomes positive). Columns (5) and (6): the 
redshift zbb and time tBB necessary to reach the Big Brake. 
The former indicates the relative size of the universe when it 
encounters the Big Brake. (The values of t» and tsB were 
computed with the Hubble parameter Ho = 73 km/s/Mpc.) 



yo 


Wo 


z* 


U 


(lOVs) 


ZBB 


t BB (lOVs) 


-0.90 


0.635 


-0.024 




0.3 


-0.068 


1.0 


-0.85 


0.845 


-0.158 




2.4 


-0.194 


3.1 


-0.85 


0.860 


-0.162 




2.4 


-0.198 


3.1 


-0.85 


0.875 


-0.166 




2.5 


-0.201 


3.2 


-0.80 


0.890 


-0.363 




6.2 


-0.390 


6.9 


-0.80 


0.905 


-0.384 




6.7 


-0.409 


7.3 


-0.80 


0.920 


-0.408 




7.2 


-0.432 


7.9 



TABLE III: As in Table HU for k = 0.4. 



yo 


wo 


z* 




(lOVs) 


zbb 


t BB (lOVs) 


-0.80 


0.710 


-0.059 




0.8 


-0.106 


1.6 


-0.80 


0.725 


-0.059 




0.8 


-0.105 


1.6 


-0.80 


0.740 


-0.060 




0.8 


-0.105 


1.6 


-0.75 


0.815 


-0.144 




2.1 


-0.184 


2.9 


-0.75 


0.830 


-0.147 




2.2 


-0.187 


3.0 


-0.75 


0.845 


-0.150 




2.2 


-0.189 


3.0 


-0.70 


0.845 


-0.241 




3.8 


-0.276 


4.6 


-0.70 


0.860 


-0.248 




4.0 


-0.282 


4.7 


-0.70 


0.875 


-0.256 




4.1 


-0.290 


4.9 


-0.70 


0.890 


-0.264 




4.2 


-0.298 


5.0 


-0.65 


0.860 


-0.358 




6.2 


-0.387 


7.0 


-0.65 


0.875 


-0.372 




6.5 


-0.400 


7.2 


-0.65 


0.890 


-0.388 




6.8 


-0.415 


7.6 


-0.65 


0.905 


-0.406 




7.2 


-0.432 


8.0 


-0.60 


0.875 


-0.521 




10 


-0.542 


11 


-0.60 


0.890 


-0.551 




11 


-0.571 


12 


-0.60 


0.905 


-0.587 




12 


-0.605 


13 


-0.55 


0.875 


-0.756 




19 


-0.766 


20 


-0.55 


0.890 


-0.837 




25 


-0.845 


26 



Finally we have evolved numerically backward in time 
some of the trajectories crossing the la domain, until 
they reached one of the Big Bang singularities of the 
model. All trajectories we have checked originate from 
the singularity at |s| = 1. In other words, they start from 
the horizontal boundaries of the rectangle in the phase 
plane (T, s), and depending on whether they evolve into 
an infinite de Sitter expansion or reach the Big Brake 



TABLE IV: As in Table [TH for k = 0.6. The evolutions into 
a Big Brake Singularity compatible with supernova observa- 
tions are more numerous with increasing k. 



yo 


Wo 


z* 


u (ioV«) 


zbb 


t BB (lOVs) 


-0.75 


0.665 


-0.039 


0.5 


-0.088 


1.4 


-0.70 


0.755 


-0.098 


1.4 


-0.145 


2.3 


-0.70 


0.770 


-0.100 


1.5 


-0.145 


2.3 


-0.70 


0.785 


-0.101 


1.5 


-0.146 


2.3 


-0.70 


0.800 


-0.102 


1.5 


-0.146 


2.3 


-0.65 


0.815 


-0.168 


2.6 


-0.209 


3.4 


-0.65 


0.830 


-0.171 


2.6 


-0.212 


3.4 


-0.65 


0.845 


-0.175 


2.7 


-0.215 


3.5 


-0.60 


0.830 


-0.240 


3.9 


-0.277 


4.7 


-0.60 


0.845 


-0.247 


4.0 


-0.283 


4.8 


-0.60 


0.860 


-0.254 


4.1 


-0.289 


4.9 


-0.60 


0.875 


-0.261 


4.2 


-0.296 


4.0 


-0.55 


0.845 


-0.325 


5.5 


-0.357 


6.3 


-0.55 


0.860 


-0.335 


5.7 


-0.366 


6.5 


-0.55 


0.875 


-0.347 


5.9 


-0.377 


6.7 


-0.55 


0.890 


-0.359 


6.2 


-0.389 


7.0 


-0.50 


0.845 


-0.411 


7.5 


-0.439 


8.3 


-0.50 


0.860 


-0.427 


7.8 


-0.453 


8.6 


-0.50 


0.875 


-0.444 


8.2 


-0.469 


9.0 


-0.50 


0.890 


-0.463 


8.6 


-0.488 


9.4 


-0.45 


0.860 


-0.533 


10 


-0.554 


11 


-0.45 


0.875 


-0.557 


11 


-0.577 


12 


-0.45 


0.890 


-0.584 


12 


-0.603 


13 


-0.45 


0.905 


-0.616 


13 


-0.633 


14 


-0.40 


0.860 


-0.658 


15 


-0.673 


16 


-0.40 


0.875 


-0.693 


16 


-0.707 


17 


-0.40 


0.890 


-0.733 


18 


-0.745 


19 


-0.40 


0.905 


-0.779 


21 


-0.789 


22 


-0.35 


0.860 


-0.814 


23 


-0.822 


24 


-0.35 


0.875 


-0.865 


28 


-0.872 


29 


-0.35 


0.890 


-0.927 


36 


-0.930 


37 


-0.30 


0.845 


-0.955 


43 


-0.957 


44 



singularity, they belong to either type II or III. 



V. CONCLUDING REMARKS 

In this paper we have shown that the tachyon cosmo- 
logical model of Ref. Q allows for a consistent set of 
trajectories which are compatible with the supernovae 
type la data. 

We have found that, among these, for positive values 
of the parameter k of the model, there is a subset of 
evolutions which end up into a Big Brake singularity and, 
for the latter, we have computed the relevant Big Brake 
parameters zbb and tBB- 
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k=-0.4 k=-0.2 




FIG. 3: (Color online) The future evolution of those universes, which are in a 68.3% confidence level fit with the supernova 
data. The la contours (black lines in the z — plane) are from Fig [5] (the parameter plane (yo,wo) is the z = plane here). 
The sequence of figures and the values of k are the same as on Fig. [2] The short and thick (blue) line in the plane of initial 
conditions separates the la parameter ranges for which the universe evolves into a de Sitter regime or towards the Big Brake 
singularity. Future evolutions towards the Big Brake singularity of the universes selected by the comparison with supernovae 
data become more frequent with increasing k. 



The compatibility of cosmological evolutions possess- 
ing soft cosmological singularities with the supernovae 
type la data was studied in Q. Curiously, it was found 
in Ref. [9] that a sudden singularity may take place in al- 



ready a very close future, even less then 10 million years. 
However this analysis was purely kinematical, and we 
also note that the parameters in our model (as given by 
the tachyonic dynamics) near the Big Brake singularity 
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fall outside the range considered in [9(. The problem of 
stability of a cosmological evolution in the vicinity of such 
singularities was studied in 10]. 

Finally, we may ask why the model proposed in [7j 
is worth studying. First, the soft (sudden) cosmological 
singularity of the Big Brake type arises in our model in 
a very natural way as a particular class of solutions of 
the dynamical system @-([7]). Second, the model has 
another interesting feature. A subtle interplay between 
geometry and matter, induces a change of the very nature 
of the latter: it transforms from tachyon into a "pseudo- 
tachyon" field (see Q for details). We point out that a 
similar effect was observed also in scalar-phantom cos- 
mological models • Phenomena of this kind represent 
a distinguishing feature of general relativity [HI : the re- 
quirement of self-consistency of Einstein equations can 
impose the form of the equations of motion for the mat- 
ter. 

Thus, in spite of being a toy model, the tachyon cosmo- 
logical model [7[ can serve as a prototype of realistic (i.e. 



compatible with observational data) cosmological models 
which may lead to a final fate of the Universe, different 
from the infinite quasi - de Sitter expansion of the ACDM 
model. What will actually happen in the future is left to 
our far away descendants to experience! 
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